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ABSTRACT 

Supernovae of Type lib contain large fractions of helium and traces of hydrogen, which 
can be observed in the early and late spectra. Estimates of the hydrogen and helium 
mass and distribution are mainly based on early-time spectroscopy and are uncertain 
since the respective lines are usually observed in absorption. Constraining the mass 
and distribution of H and He is important to gain insight into the progenitor systems 
of these SNe. 

We implement a NLTE treatment of hydrogen and helium in a three-dimensional 
nebular code. Ionisation, recombination, (non-)thermal electron excitation and Ha 
line scattering are taken into account to compute the formation of Ha, which is by 
far the strongest H line observed in the nebular spectra of SNe lib. Other lines of 
H and He are also computed but are rarely identified in the nebular phase. Nebular 
models are computed for the Type lib SNe 1993J, 2001ig, 2003bg and 2008ax as well 
as for SN 2007Y, which shows Ha absorption features at early times and strong Ha 
emission in its late phase, but has been classified as a SN lb. We suggest to classify 
SN 2007Y as a SN lib. Optical spectra exist for all SNe of our sample, and there is 
one IR nebular observation of SN 2008ax, which allows an exploration of its helium 
mass and distribution. 

We develop a three-dimensional model for SN 2008ax. We obtain estimates for 
the total mass and kinetic energy in good agreement with the results from light- 
curve modelling found in the literature. We further derive abundances of He, C, O, 
Ca and 56 Ni. Estimates of the H mass are difficult but some constraints are derived. 
We demonstrate that Ha absorption is probably responsible for the double-peaked 
profile of the [O i] A A 6300, 6363 doublet in several SNe lib and present a mechanism 
alternative to shock interaction for generating late-time Ha emission of SNe lib. 
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INTRODUCTION 

Massive stars (> 8 M©) collapse when the nuclear fuel in 
their central regions is consumed, producing a core-collapse 
supernova (CC-SN) and forming a black hole or a neutron 
star. CC-SNe with a H-rich spectrum are classified as Type 
II. If the envelope was stripped to some degree prior to the 
explosion, the SNe are classified as Type lib (strong He lines, 
and weak but clear H), Type lb (strong He lines but no H), 
and Type Ic (no He or H lines) (|Filippenkoll 19971 ). 

Large amounts of hydrogen (^> 1 M0) cause the pro- 
genitor to become a red super-giant prior to explosion, 
which in turn leads t o a broad peaked SN light curve (e.g. 
iGrassberg et al.lll971f ). Since this is not observed in SNe of 
Type lib there is an upper limit for the hydrogen mass of 



less than one solar m ass (|Nomoto et al.lll993 ; IUtrobinlll994l ; 
IWooslev et al.l l 1994). It is not clear how SNe lib manage to 
keep just a thin layer of hydrogen. 

While massive (and therefore hot, O-Type) stars (M > 
25 Mq) can blow o ff their hydrogen envelo pe by radiatively 
driven winds (e.g. lEldridge fe Vinkl I2006T) . mass estimates 
for the progenitors of SNe lib (e.g. IWooslev et al.l 1994; 
Mazzali etall 120091 : ISilverman et al.l 12009) ; lHamuv et all 



20091 ) suggest that they may not be massive enough to lose 



most of the H envelope via this process. Binary interaction 
can cause mass transfer between two interacting stars and 
would allow a Type lib progen itor to lose most of its hy- 
drogen (e.g. lWooslev et al.lll994J ). There is o bservational ev- 
idence that the progenitors of SNe 1993J jIAldering et al.l 
1 1994 iMaund et all 12004 iMaund fc Smarttl 120091 ), 2001ig 



© 0000 RAS 



Maurer et al. 



l|Rvder et alj|2006h [also see lKotak fc Vinkl l|2006l ) for an al- 
ternative interpretati on of the SN 2001ig data] and 2008ax 
l|Crockett et al.ll2008l ) may have been part of binary systems. 
The re is evidence of shock interaction from nebular 
Haje.g. IChugailll99ll : IPatat et al.lll995l : iHouck fc Franssonl 
1 19961 . also see this wor k) as well as from radio 



e.g. Fransson fc Biornssonl 1 19981 ; ISoderberg et al.l 120061 ; 



Chevalier fc Soderberd 20101) and X-ray observations (e 



Chevaliej|l98ll; ISoderberg et alJlJOOrJ ; iNvmark et al.|[200£ 
Chevalier fc Soderberdl2010r ). suggesting that at least som 



SNe lib are surrounded by massive stellar winds. This wind 
could also be produced by the massive companion star in 
the binary scenario. 

At early phases of SNe lib the energy emitted by 
Ha is probably provided by the radioactive decay chain 
56 Ni -> 56 Co -> 56 Fe. However, it has been show n for 
SN 1993J (|Patat et al.lll995l ; IHouck fc Franssonlll996l ) that 
~ 150 days after explosion this mechanism becomes inef- 
fective suggesting that an additional source of energy is re- 
quired to explain the observed Ha luminosities (however this 
may not be true; se e Section [5]) . A similar result had been 
obtained for SNe II (|Chugailll99ll ). This is further confirmed 
by the detec tion of a flattening of t he Ha luminosity decay at 
late phases (JMatheson et al.l l2000). This could be inconsis- 
tent with radioactive decay, but may be explained by shock 
interaction. 

Shock interaction can also be detected through ra- 
dio and X-ray observations. As the SN ejecta propa- 
gate into the circumstellar medium, they are deceler- 
ated, creating internal energy which is dissipated by ra- 
diative processes. There are several solu tions for differ- 
ent scenarios of shock interaction (e.g. [C hevalier Il98ll ; 
ISuzuki fc Nomotol Il995l ; iFransson fc Biornssonl Il998f l. but 
a clear interpretation of the observations is often difficult. 
Micro-physical processes (e.g. the formation of magnetic 
fields) are poorly understood, and it is not clear in de- 
tail how the energy released by the shock is transferred 
into radiation. The SN envelope and circumstellar den- 
sity profiles are not known and have to be treated as free 
parameters. Several authors found evidence for inhomoge- 
neous wind structures and deviations fro m wind-like den- 
sity profiles of the e xternal medium (e . g. [ Franssonl 1 1994 ; 
Ivan Dvk et ail Il994l ; ISuzuki fc Nomotol 1 19951 ). This may 
however be an artefact of an inaccurate tre atment of shock 
physics l|Fransson fc Biornssonl 1 19981 , 120051 ). A very careful 
treatment of shoc k interaction seem s necessary to obtain re- 
liable results fe.g. IFransson fc Biornssonl 1 19981 ). 

Asphericities in the inner and outer ejecta are evi- 
dent in at least some CC-SNe. Three indicators are ve- 
locity diffe rences of Fe and l ighter element lines at late 
time s fe.g. IMazzali et al.ll200ll ). polarisation measurements 
(e.g. lHonicblll99rl ) a nd double-peaked or asymmetric emis- 

line profiles (e.g. IMazzali et al.ll2005l;lMaeda et alf2 008; 
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Modiaz et"aLll2008l;lTaubenberger et alJ|2009l ; IMaurer et al.1 

2010) [also sec Milisavli evic et al.l (2010) for an alternative 
interp retation, which is however doubtful (JMaurer et al.l 
|2010| )]. Indirect evidence also emerges from a compariso n 
of the inner and outer ejecta velocities (jMaurer et al.ll2010h . 
We find evidence that SN 2008ax is a non-spherical event 
(see Section [3} . 

Between 100 and 200 days after explosion most SNe 
enter their nebular phase. Owing to the decrease of den- 



sity the SN ejecta become transparent at most wavelengths, 
which makes it possible to observe the innermost parts of 
the SN. Nebular modelling provides the opportunity to de- 
rive information about the SN ejecta mass, kinetic energy, 
abundances and geometry, which can hardly be obtained by 
other methods at least for the central regions. 

Nebular observations are available for very few SNe lib. 
To our knowledge there are nebular spectra of SNe 1993J, 
2001ig, 2003bg, 2006T, 2008aq and 2008ax. The quality of 
the SN 2006T spectrum is poor. For SN 2008aq there is no 
light curve and therefore there are no flux-calibrated spectra. 
Therefore, neither SNe are included in our analysis. 

Although early time Ha absorption features are 
observed, SN 2007Y was classified as a SN lb by 
IStritzinger et al.l 1)20091 ). who argued there is very little 
hydrogen in the envelope and that the late Ha emis- 
sion is powered by shock int eraction. On the other hand, 
IChevalier fc Soderberd (J2010J) found that the Ha emission of 
SN 2007Y must be powered by radioactive decay before 300 
days after explosion. SN 2007Y is included in our analysis to 
investigate this contradiction and a possible re-classification. 

In Section 2 we describe the implementation of hydro- 
gen and helium in our nebular code. In Section 3 we de- 
scribe a three-dimensional model of SN 2008ax, give some 
estimates for the total mass and kinetic energy and compare 
them to results available in the literature. In Section 4 we 
present nebular models for the SNe 1993J, 2001ig, 2003bg, 
2007Y. Since nebular models for these SNe are available in 
the literature already and since we cannot determine their 
helium density we concentrate on their Ha emission. In Sec- 
tion [5] we present a scenario alternative to shock interaction 
for explaining strong late-phase Ha emission. In Section [6] 
our results are discussed. 



2 H AND HE IN THE NEBULAR PHASE 

In this section we describe the im plementation of hydro- 
gen and helium in our nebular code dRuiz-Lapuente fc Lucvl 
I1992I ; IMazzali et ail l200ll , I2007J ; IMaurer et al.ll2010l ). The 
code computes the energy deposition from the radioactive 
decay of 56 Ni and 56 Co in the SN ejecta, uses the energy de- 
posited to compute ionisation and excitation, and balances 
this with gas cooling via line emission in a NLTE sc heme. 
The method follows that developed bv lAxelrodl (|l980h . 



2.1 Hydrogen 

A full NLTE treatment of hydrogen is implemented in our 
three-dimensional nebular code. We obtain radiative transi- 
tion and ionisation rates for the hydrogen at om from analyt - 
ical solutions available in the literature fe.g. lBurgesslll965l ). 
Collisional rates for allowed transitions between excited lev- 
els n > 3 are uni mportant, but a re included using an approx- 
imation (|van Regemorterl 1 19621 ) . Collisional r ates between 
quant u m levels n = 1,2 and 3 are taken from IScholz et al.l 
(| 199ah : ICallawavl (I1994T ). The 2s level is additionally con- 
nected to the ground state by two-photon decay, and the 2s 
and 2p levels are coupled by electron and proton collisions. 
Hydrogen is ionised by non-thermal electrons (pro- 
duced by Compton scattering of 56 Co radioactive decay 7- 
radiation) and by UV-radiation emitted by helium if hy- 
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drogen and helium are mixed. The electron impact ionisa- 
tion rate can be obtained comparing the atomic and elec- 
tronic loss-functions with the ionisation cross-sections (e.g. 
lAxelrodl Il980l ; iMaurer fc Mazzalil l201Ch . Thermal electron 
collisional ionisation and recombination is not important. 
Excitation by non-thermal electrons is important. We com- 
pute non-t hermal electron excita tion rates using the approx- 
imation of iRozsnvai et al.l (|19801 ). 

We compare our hydrogen and helium ionisation and 
non-thermal excitation rates to calculations from Hachinger 
et. al. (in prep.) who calculate non-thermal ionisation and 
excitation rates solving an energ y balance equation derived 
from the Spencer- Fano equation (|Xu fc McCravlll99ll ; ILucvi 
ll991T ). Although these methods are quite different, the ther- 
mal excitation and ionisation rates at various electron and 
atomic densities agree to 10% — 20% for both hydrogen and 
helium. 

For levels with principal quantum number n > 5 we 
only consider recombination (these levels are almost com- 
pletely depopulated and have no influence on the nebular 
spectrum), while a full NLTE (de-)excitation treatment is 
performed for lower levels. 

The n = 2 level of H can effectively scatter the back- 
ground r adiation field. This is a v ery important process in 
SNe lib |Houck fc Franssonl Il996h [Hou96] since Ha comes 
into resonance with the [O i] A A 6300, 6363 doublet, which 
carries an important fraction of the nebular flux. The optical 
depth of transitions between levels njs 3 is too low to cause 
observable line scattering at late epochs. 

Clumping of hydrogen is usually neglec t ed in SNe lib 
i Patat et all 1 19951 ; iHouck fc Franssonl 1 19961 ; iMazzali et all 



12009J), but it can increase the optical depth of Hq consider- 
ably. Under certain conditions, as we show in Section [5] it 
can also increase the Ha emission significantly. 



2.2 Helium 

Accurate atomic data for He are available in the liter- 
ature (e.g. see TIPTOPbase 0, NIST 0). For technical 
reasons it was more convenient for us to calculate the 
He radiative transit ion and ionisation rates using a quan- 
tum defect method jBates fc Dam gaard 19 491; ISeatonll 19581 ; 



3fT95l; 

i (2010) 



iBurgess fc Seatonl (19601 ') [also see IMaurer fc Mazzalil (|201C 
for an application to oxygen] , which is very accurate for he- 
lium. The disagreement between the calculations and NIST 
recommended data is a few percent in the worst case, having 
no observable influence on the line calculations. 

Thermal electron excitation rates for all quantum levels 
n < 4 are taken from iBerrington fc Kingston! |l_987). The 
2s( 1 S) level is conne cted to the ground state by two-photon 
emission (2PE) (e.g. iDrake et al.lll969T ). which turns out to 
be important in the nebular phase of SNe of Type lib (in 
the early phase radiative excitation of the 2p levels allows 
effective de-excitation via the 2p — > Is transition, which can 
be more important than 2PE then). 

Helium is mainly ionised by non-thermal electrons (see 
hydrogen). The recombination rates into all levels n ^ 4 are 
calculated by simulating recombination into and cascading 



1 http://cdsweb.u-strasbg.fr/topbasc/ 

2 http://www.nist.gov/index.html 
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Figure 1. Illustration of the three-dimensional model of SN 
2008ax, described in the text. Regions drawn in red colour are 
rich in heavy elements and contain no helium. Regions indicated 
in blue are helium-rich. The outermost region, indicated in green, 
is hydrogen-dominated. The x-axis is associated with the equator. 
In the centre (below 4500 km s — 1 ), heavy elements (red) are con- 
centrated towards the pole, while ejecta containing helium (blue) 
are concentrated in the equatorial regions. 



from excited levels n ^ 30. Levels with n ^ 4 are treated 
in full NLTE. Non-thermal electron excitati on is important 
and is included using the approximation of IRozsnvai et al.l 
l|l980h . 

Helium can be observed in emission during the nebular 
phase of SNe lib since its mass is large and since it is dis- 
tributed down to a few thousand km s _1 . Owing to the large 
excitation potential of He (~ 20 eV), an important fraction 
of the deposited energy can be radiated away by other ele- 
ments which may be present in small fractions. This means 
that the electron temperature decreases and the He lines be- 
come weaker while the other elements emit more strongly. 
Since all optical He lines are blended with lines only the IR 
He 1 10830 A and 20587 A lines can be identified clearly 
(at least in SN 2008ax) , which allows an estimate of the He 
mass and distribution. Since both lines are optically thick 
around 100 days after explosion, they can be influenced by 
line scattering. 



3 SN 2008AX 

Nebular spectra of SN 2008ax are available at 131, 149, 
266, 280 and 359 days after explosion. The spectrum at 131 
days covers the range between 9000 and 25000 A. The other 
spectra cover a range between 5000 and 10000 A. For the 
modelling we use a distance modulus of fi — 29.92 ± 0.29 
mag and an extinction of E(B — V) — 0.4 ± 0.1 mag (see 
iTaubenberger et al.|[2jnjj)[Tau l0l. Also see iPastorello et all 
(20081; IChornock et alj (|2010T ) for observations and discus- 
sion of SN 2008ax. 



© 0000 RAS, MNRAS 000, 000-000 



Maurer et al. 



Table 1. SN 2008ax. (A) Best fit three-dimensional model (B) All 56 Ni confined below 2800 km s" 1 ; 
inconsistent with nebular Fe-group line observations. Hydrogen is not included in the total mass and kinetic 
energy estimate in both models but should be of the order of 0.1 Mq and 10 50 ergs. 
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Table 2. Properties of SN 2008ax obtained by different methods: (A) Multi-dimensional nebular modelling 
(B) one-zone light-curve calculation (C) semi-analytical light-curve modelling (D) numerical light-curve mod- 
elling using the radiation transport code STELLA on a SN lib explosion model (E) light-curve comparison 
of SN 2008ax and SN 1993J. 

*Please note that (C), (D) and (E) are computed for E(B — V) = 0.3 mag instead of 0.4 mag. In rough 
approximation a comparison can be made by multiplying all quantities of (C), (D) and (E) by a factor of 
1.15. These values are given in brackets. 
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There is observational evidence for asymmetry in SN 
2008ax. The double-peaked profiles of the He I IR 10830 and 
20587 A lines can be i nterpreted as a t o rus-sh aped distribu- 
tion of helium [TaulO lChornock et al.l (|2010T )] . In addition, 
the blue wing of these lines is stronger, which can be inter- 
preted as an asymmetry along the line of sight, but may also 
be a scattering effect. We place the observer in the equato- 
rial plane, which is necessary to produce the double-peaked 
He profiles in this model. 

The choice of this geometry (see Figure [l] for illustra- 
tion) is motivated by the observations but is probably not 
unique. Different geometries may reproduce acceptable fits 
to the observations as well. 

In contrast to the He IR lines (see Figure [2j, most 
other lines are single peaked (see Figure[3]). This is expected, 
since heavy elements are concentrated in the core. An ex- 
ception is the [O i] AA 6300, 6363 doublet. While we agree 
that the profile of these lines may be shaped by geometr 
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[TaulO] in other types of CC-SNe (e. 
Maeda et a.l.ll2008J: [Modiaz et alJlJOO 
20091 ; iMaurer et alj|201fj ), we think that in SNe lib Hq line 



I 



Taubenberger et al.l 



scattering is responsible for the splitting of the [O i] AA 6300, 
6363 doublet (also see Section [4] for other SNe lib). 

The velocity of the Ha absorption minimum saturates 
at ~ 12500 km s _1 about 40 days after explosion. This 
is usually interpreted as the lower boundary of hydrogen. 
Since hydrogen in lower layers may be strongly ionised (see 
Appendix A), this estimate is somewhat uncertain. It may 
rather be the lower boundary of H I but not of H II. 

We can reproduce the double-peaked profile of the [O i] 
A A 6300, 6363 doublet very well placing less than 0.1 M of 
hydrogen between 12000 and 12500 km s _1 (see Figure [3]). 
An exact estimate of the H mass is not possible since shock 
interaction or clumping may be important in this region. 



This can explain why other lines, like [O i] 5577 A are single- 
peaked. This is discussed in more detail in Section [6] 

Our model consists of 128 angular and 8 radial cells 
(spherical geometry; 1024 cells in total). The radial cells 
have outer boundaries at 2000, 2800, 4500, 6000, 6600, 9500, 
12000 and 12500 km s _1 . This choice relates to the physical 
properties of our model for SN 2008ax. 

Below 2000 km s _1 the ejecta are spherically symmet- 
ric. This innermost part of the SN is dominated by 56 Ni, 
oxygen and calcium. The total mass and mass fractions in 
this zone are quite uncertain, but their contribution to the 
global properties of the SN is small. No He is present (< 
10" 4 M ). 

The region between 2000 and 2800 km s _1 is similar 
to the innermost part, but some cells, preferentially in the 
equatorial plane, contain small fractions of He (~ 10 -3 M ). 
The bulk of the oxygen, calcium and carbon of SN 2008ax 
is located in this region. In contrast to calcium, carbon is 
probably not present below 2000 km s _1 . The confinement of 
C to a thin low v elocity layer is in agr eement with theoretical 
predictions (e.g. iNomoto et al.lll993l ). 

The region between 2800 and 4500 km s _1 is domi- 
nated by oxygen and 56 Ni in some cells and by He in others. 
This separation improves the reproduction of the observa- 
tions. On the one hand He emission can be observed down 
to 2000 km s _1 and the He IR line ratio suggests that there 
is no strong mixing with heavy elements (which decreases 
the electron temperature and influences the strength and 
ratio of the He I 10830 and 20587 A lines). On the other 
hand, the iron lines observed are much broader than 2800 



km 



Placing all the Ni mass necessary to reproduce 



the observed flux below 2800 km s _1 causes too narrow and 
too strong Fe emission lines. A good compromise is found by 
allowing a separation of the 56 Ni and He-rich ejecta. Large- 
scale structures of burned and unburned material are not 



© 0000 RAS, MNRAS 000, 000-000 



Hydrogen and helium in SNe lib 5 



Velocity [10 km s 1 
+10 -10 




1.03 1.08 1.13 

Wavelength \fj.] 



1.95 2.05 

Wavelenqth [/x] 



Figure 2. IR spectrum of SN 2008ax at 131 days after explo- 
sion. The He I 10830 A line is shown on the left and the He I 
20587 A line is shown on the right. The observations are shown 
in black while the synthetic flux is shown by the dashed red line. 
The synthetic flux is obtained using the three-dimensional model 
described in the text. Since almost no background flux is pro- 
duced by our nebular code near 20000 A we added some constant 
flux at the level of the observations between 19000 and 22000 A 
to allow He I 20587 A line scattering. In contrast to the optical 
He lines, the 10830 A and 20587 A lines can be identified clearly, 
which allows a determination of the He density field. Two kinds 
of asymmetry can be observed. First, the double-peaked nature 
suggests a torus-shaped distribution of He. Second, the blue side 
of the profiles of both lines is stronger, which suggests an asym- 
metry along the line of sight. 



M 










25 






i 


- 


20 










- 


15 








I Ir I 


- 


10 










- 


5 


\/ v>SJ '\ Ml J 


A 


r 


\ J I / \ 


A 


n 






v-/ *. ' \ „ r . * _| 


J\^. 



5000 6000 7000 8000 

Wavelength [A] 



Ha Velocity [10 km s M 
15 +10 +5 



9000 



10000 




6300 6600 

Wavelength, [Al 



6900 



unexpected depending on the explosion scenario. The He 
rich cells contain some fraction of carbon, oxygen, sodium 
and calcium (in total ~ 30%). These elements are strongly 
excited by the energy absorbed by the He layer. 

The material above 4500 km s _1 is dominated by he- 
lium, which constitutes most of the mass of SN 2008ax. The 
He rich layer reaches out to at least ~ 9500 km s -1 . The 
composition and density of the zone between 9500 and 12000 
km s _1 is not clear. The absorption profile of the He I IR 
lines suggests that there is no He I 20587 A line scattering 
above ~ 10000 km s _1 (see Figure [2J. An upper limit to 
the He mass in this region is ~ 0.05 Mq. The observation 
of the [O i] AA 6300, 6363 doublet shows that there is no 
strong Ha line scattering below 12000 km s _1 (see Figure 
[3]). A possible solution may be that there is some He (of 
the order of 0.1 Mq) in this region, but mixed with small 
fractions of neutral elements, like H I. In this case contin- 
uum destruction of He I 584 A photons could red uce the 
optical depth of the He I 20587 A line significantly (|Chugail 
1 19871 : iLi fc McCravl Il995l ) . This effect strongly depends on 
the composition of the He layer. However, since the mass in 
this region is expected to be low, the effect on our total mass 
estimate would presumably be less than 10%. The estimate 
of the kinetic energy could be influenced more strongly. 

The 266, 280 and 359 day Ha observations cannot be 
reproduced using any reasonable amount of hydrogen (of 
the order of 0.1 Mq). It seems likely that some additional 
mechanism of Ha emission becomes important between 150 
and 200 days after explosion, as in SN 1993 J [Hou96]. In 



Figure 3. Optical spectrum of SN 2008ax at 149 days after ex- 
plosion. Observations are shown in black while the synthetic flux 
is shown by the dashed red line. The synthetic flux is obtained us- 
ing the three-dimensional model described in the text. Although 
the He I IR lines are shaped by the torus-like He distribution, all 
other lines are single peaked, since the heavy elements are con- 
centrated in the core. The (double) peak of the [O i] AA 6300, 
6363 doublet profile is caused by He* absorption between ~ 12000 
and 12500 km s — 1 . The Ho absorption minimum is indicated by 
the vertical dashed line. There seems to be continuum flux, espe- 
cially between 7000 and 10000 A, which is not reproduced. This 
is because of the epoch (149 days), which is too early for a strictly 
nebular treatment. 



Section [5] we propose a mechanism, based on mixing and 
clumping of hydrogen and helium. Shock interaction is also 
a possibility and is discussed in Section [6] 

The O I 7774 A line is reproduced well at 266 days and 
later bu t is too weak at 149 days after explosion, which is ex- 
pected IjMaurer fc MazzalilbOlOT ) since there is no clumping 
of the ejecta in our model. 

The He I 10830 A line may contain some contribution 
from Si and S lines. However, in SN 2008ax this contribu- 
tion is weak. It seems likely that He I 20583 A and possibly 
He I 10830 A are influenced by line scattering. Both lines 
have optical depths larger than one around 130 days after 
explosion. However, it is not clear if there is enough flux 
that can be scattered. Our nebular code produces almost no 
emission on the blue side of both lines. However, at least for 
He I 20587 A this is in conflict with the observations, which 
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show some background flux around 20000 A. Since the spec- 
trum at 131 days is probably not completely nebular, it is 
not surprising that certain features of the spectrum are not 
reproduced well (also see the background flux in the optical 
spectrum at 149 days). 

To handle this problem we introduce some artificial 
background flux between 19000 and 22000 A at the flux level 
of the observations. This allows the He I 20587 A line to in- 
crease by line scattering. However, since this background is 
completely artificial, one can not expect that the modelling 
is accurate. 

The masses estimated for He, C, O, Ca and 56 Ni are 
shown in Table 1. The uncertainties on these estimates 
should be of the order of (several) 10%. Other elements 
carry larger uncertainties. We also compute the total mass 
and kinetic energy, excluding hydrogen. Our best-fit three- 
dimensional model is listed Table 1, row (A). Since some 
SN lib m odels predict a confi nement of 56 Ni to quite low 
velocities (No moto et al.ll 19931 ). we additionally show an es- 
timate obtained by placing all 56 Ni below 2800 km s _1 (B), 
which is, however, in conflict with the Fe-group emission line 
observations. 

A comparison of the total mass of models (A) and 
(B) shows that the uncertainties owing to the distribution 
of 56 Ni is about half a solar mass. However, since model 
(B) stands in clear contradiction to the nebular Fe-group 
emission line observations model (A) has to be preferred. 
There are also uncertainties from He I background scatter- 
ing, clumping and from the atomic data. Additionally, we 
may underestimate the mass in the outer regions. 

Owing to all the uncertainties described above, we esti- 
mate a total mass between 2.2 and 3.2 Mq and a 56 Ni mass 
between 0.10 and 0.12 Mq. A lower limit on the kinetic en- 
ergy is 0.9 TO 51 ergs (for a total mass of 2.7 Mq). Since 
small amounts of high velocity H and He could increase the 
total kinetic energy considerably we estimate a total kinetic 
energy between 0.7 and 1.7 • 10 51 ergs (for a total mass 
between 2.2 and 3.2 M ). 

In addition there is some uncertainty (~ 40%) owing 
to extinction (we used E(B — V) — 0.4 ±0.1 mag) and 
distance (9.6 ± 0.3 Mpc). Since these uncertainties would 
influence the total mass and the 56 Ni mass estimate simul- 
taneously we assume that the uncertainty owing to distance 
and extinction is ~ 20% on each of these quantities. Thus 
we estimate (including 0.1 Mq of hydrogen with a kinetic 



energy of 10 ergs) a total mass of 2.i 



M Q , a total 5b Ni 



and a total kinetic energy l.O^o'g ' 10 51 



mass of O.lOl^;^ 
ergs. 

In Table 2 we compare our estimates to results from 
other groups. [TaulO] obtain estimates from one-zone- 
modelling of the light curve. They estimate a total mass 
of ~ 4.6 Mq, a kinetic energy of 6 ■ 10 51 ergs and a 56 Ni 
mass of 0.1 Mq. The estimates of the total mass and kinetic 
energy are larger than found in this work. The uncertain- 
ties estimated by [TaulO] are large and the re sults agree 
within these uncertainties. iRoming et alj (|2009r ) model the 
light-curve of SN 2008ax using a combination of an analyti- 
cal light-curve model and a Monte-Carlo routine. Our results 
ar e roughly consi s tent, however the kinetic energy estimated 
bv lRoming et al . (2009) seems to be too low in general. This 
has also been found by [TaulO]. 

A comparison to lTsvetkov et al.l 1)20091 ) shows that our 



results agree rather well. Since the extinction w as estimated 
to be lower in the work of lTsvetkov et al.l (120091 ) we multiply 
their results by a factor of 1.15 for comparison. This gives 
0.127 M Q of 56 Ni, a total mass of 2.7 M and a kinetic 
energy of 1.7- 10 ergs. We may underestimate the kinetic 
energy but our results agree within the errors. 

Our estimates for the 56 Ni and total mass, as well as 
for the kinetic energy are also consistent with the results of 
iPastorello et all (|2008D . 



4 OTHER SNE OF TYPE IIB 



Nebular model s 
(Silverman et al. 



for SNe 199 3J [Hou96], 2001ig 
20091). 2003bg (JMazzali et all 12009b 



and 2007Y ( Stritzinger et al.l 120091 ) exist but a treatment 
of helium was not possible, since there are no IR nebular 
spectra ([Hou96] included He in their analysis but had 
no IR observations for comparison). An exact treatment 
of hydrogen is also difficult since it is observed as Ha, 
which may be influenced by scattering, clumping and shock 
interaction. 

Therefore, we restrict our analysis to the question: can 
Ha be powered by radioactive energy deposition, as de- 
scribed by [Hou96]? The answer to this question is less 
model dependent than estimating element masses or distri- 
butions. At late epochs (> 200 days) the SN ejecta are illu- 
minated homogeneously by 7-radiation (since the 7-optical 
depth is low) and asymmetries in the 56 Ni or hydrogen dis- 
tribution are not very important. In addition we can ob- 
tain some constraints on the H and He distribution from 
early and late-time line-width observations. The maximum 
H mass is restricted by light c urve observations and should 
be less than one solar m ass (|Nomoto et alj 1 19931 ; lUtrobinl 
ll994l ; IWooslev et al.lfl994T ). 

Early time observations show that all SNe of our sample 
have absorption minima of Ha between 13000 and 20000 
km s _1 during the first days after explosion. The velocity of 
these absorption minima decreases rapidly in the following 
tens of days and bottoms out at ~ 10000 — 13000 km s _1 in 
all SNe around 30 — 40 days after explosion. This saturation 
behaviour is usually interpreted as the lower boundary of 
the hydrogen layer. However, there may be some hydrogen 
at lower velocities if it is ionised completely in the early 
phase of the SN (see Appendix A). We show that there is 
an intriguing link between the velocity of the early time Ha 
absorption minimum and the profile of the [O 1] AA 6300, 
6363 doublet for all SNe of our sample (see Table 3). 

Whether a double-peaked [O 1] AA 6300, 6363 doublet is 
formed or not, depends on the velocity of the optically thick 
Ha . The optimum Ha velocity to produce a double-peaked 
O [1] profile is ~ 12000 km s _1 since then the absorption 
minimum of Ha for 6300 A lies right at the centre of the [O 
1] A 6300 line profile, which corresponds to the centre of the 
SN (bulk of oxygen). 

4.1 SN 1993J 

Nebular spectra of SN 1993J have been modelled before, 
investigating the formation of the Ha line in detail [Hou96] . 
We repeat this analysis to test whether our results are in 
agreement with previous findings. We have spectra at 118, 
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Table 3. The velocity of the absorption minimum of Ha at ~ 40 
days after the explosion (i> p hot,Ha) taken from the literature and 
the velocity of the Ha derived from the [O i] AA 6300, 6363 doublet 
profile in this work (f ne b Ha)- A physical connection between both 
velocities seems likely, confirming that the profile of the [O i] AA 
6300, 6363 doublet is strongly influenced by Ho absorption. 



Ha Velocity [10* km 
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Figure 4. Optical spectrum of SN 1993J at 206 days after explo- 
sion (black line). The synthetic flux (red dashed line) is produced 
using 0.2 Mq of hydrogen distributed between 7000 and 10000 
km s — 1 . Ha is not reproduced with sufficient strength. The [O i] 
AA 6300, 6363 doublet is single peaked, since hydrogen is located 
below 11000 km s" 1 . 



Figure 5. Left: Optical spectrum of SN 2001ig at 340 days af- 
ter explosion. Right: Optical spectrum of SN 2003bg at 301 days 
after explosion. The band shown (6100 to 6900 A) is dominated 
by the [O i] AA 6300, 6363 doublet. The observations are shown 
by the black lines. The red dashed lines show the synthetic flux 
obtained using ~ 0.3 Mq of hydrogen and a one-dimensional 
model. Hydrogen causes some weak line scattering of the oxygen- 
dominated flux around 6300 A, creating the double peaked profile 
of the oxygen doublet, but does not provide enough Ha flux to ex- 
plain the observations around 6560 A. The double peaked oxygen 
profile is not caused by geometry, but by Ha absorption around 
13200 km s- 1 (SN 2001ig) and 12800 km s^ 1 (SN 2003bg). The 
Ha absorption minima are indicated by the vertical dashed lines. 



a shock interaction dom inated phase (e.g. IPatat et al.lll99a : 
iHouck fc Fransson|[l996l ). 

The [O i] AA 6300, 6363 line in SN 1993J is not double- 
peaked. As we have shown for SN 2008ax, the double-peaked 
oxygen profile is possibly the result of Ha absorption. Since 
the optically thick Ha of SN 1993 J has a velocity < 11000 
km s _1 (see Table 3), the oxygen line is not split by the 
Ha absorption minimum. This is discussed in more detail in 
Section [6] 



172, 206, 237, 256, 300 and 363 days after explosion. While 
the 118 day spectrum is not strictly nebular, all the later 
ones are. The spectra cover a range between 4000 and 10000 
A but no IR observations are available. For the modelling we 
use a distance modulus of fi = 27.72 mag and an extinction 
of E(B -V) = 0.18 mag [Hou96]. 

[Hou96; see their Figure 1] have shown that a model 
consisting of a heavy element core (0 — 3400 km s _1 ), a 
He layer (3400 — 7800 km s _1 ) and a hydrogen dominated 
layer (> 7800 km s" 1 ) containing ~ 25% of He can produce 
synthetic spectra consistent with the nebular observations 
of SN 1993 J (except Ha at late epochs). 

Using a model similar to the one presented by [Hou96] 
we can reproduce the evolution of the heavy element lines of 
SN 1993 J at all epochs between 118 and 363 days (e.g. see 
Figure |4]| . T he hydrogen mass of SN 1993J was estimated t o 
be 0.2 M (jWooslev et al.lll994l : lllouck fc Franssonlll996h . 
Using this hydrogen mass we come to the same conclusions 
as [Hou96]. The synthetic Ha flux and the observations be- 
come inconsistent between 150 and 200 days after explosion. 
This is usually interpreted as the time of the transition to 



4.2 SNe 2001ig & 2003bg 

For SN 2003bg there is a H mass estimate of ~ 0.05 Mq ob- 
taine d from early-time absorption modelling IjMazzali et al.l 
|2009j), while there is no estimate of the H mass of SN 2001ig. 

We have nebular spectra of SN 2001ig at 256, 309 and 
340 days after explosion. The spectra cover a range between 
4000 and 10000 A. For the modelling we use a distance mod- 
ulus of fi = 30.5 mag and an extinction of E(B — V) — 0.011 
mag (|Silverman et al.ll2009T ). 

For SN 2003bg we have nebular spectra at 264 and 301 
days after explosion. The spectra cover a range between 4000 
and 10000 A. For the modelling we use a distance modulus 
of n = 31.68 mag an d an extinction of E(B — V) = 0.02 mag 
JMazzali et al.ll2009l ). 

Although the nebular Ha flux is weak in both SNe, we 
find that it cannot be reproduced using reasonable amounts 
of hydrogen. In contrast to the Ha emission, the absorption 
can be reproduced well. The [O 1] AA 6300, 6363 doublet 
is often used to investigate asymmetries of CC-SNe cores. 
However, in the case of SNe lib there is strong evidence 
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Figure 6. Optical spectra of SNe 1993J (256 days; red line), 
2007Y (248 days; black line) and 2008ax (266 days; green line). 
The spectra are scaled by arbitrary constants. Apart from the 
[O i] AA 6300, 6363 doublet and the [Ca II] emission the spec- 
tra agree extremely well. The ratio of Ha to oxygen emission 
is strongest in SN 2007Y, which is surprising since this SN was 
classified as SN lb, while SNe 1993J and 2008ax were classified 
as Type lib. SN 2007Y has the strongest nebular ratio of Ha to 
oxygen (and total) flux ever detected in a stripped CC-SN. The 
[O i] AA 6300, 6363 doublet is single-peaked, which is expected if 
the bulk of hydrogen is located below 11000 km s~ 1 . 



that the profile of the oxygen line is not shaped by geometry 
alone but also by absorption (see Section 3). To reproduce 
the oxygen profiles of SN 2001ig and 2003bg we have to place 
~ 0.3 M of hydrogen between 13000 and 13500 km s _1 (SN 
2001ig) and 12600 and 13000 km s _1 (SN 2003bg) which is 
roughly consistent with the early-time absorption minimum 
of Ha (see Table 3) of these SNe. 

Although we can reproduce the absorption features well, 
it seems unlikely that 0.3 Mq of hydrogen are necessary 
to reproduce the observations if there is some additional 
mechanism of Ha excitation (e.g. shock interaction). This 
may be expected, since Ha emission is underestimated in 
our models. 



4.3 SNe 2007Y 

Although Ha absorption was detec ted in its early spectra , 
SN 2007Y was classified as a SN lb (|Stritzinger et alJl200Sh . 
While these Ha features would be sufficient to classify SN 
2007Y as a Type lib, there is more evidence for the presence 
of hydrog en since there is s t rong Ha emission in the nebu- 
lar phase. IStritzinger et al.l (|2009l ) argue that the emission 
is cau sed by shock interaction, but IChevalier fc Soderbera 
IJ2010TI claim that the circumstellar density of SN 2007Y is 
too weak to produce the Ha luminosities observed before 
300 days after explosion. 

The nebular spectra cover a range between 4000 and 
10000 A. For the modelling we use a distance modulus of /i 
= 31.43 mag and an e xtinction of E(B — V) — 0.112 mag 
(jStritzinger et al. 2009). There is no IR nebular spectrum. 

The Ha flux can be reproduced neither at 248 nor at 288 
days using reasonable amounts of hydrogen. We therefore 
conclude that the Ha flux of SN 2007Y cannot result from 



Figure 7. Optical spectrum of SN 1993J at 363 days after ex- 
plosion (black line) . The synthetic flux (red dashed line) was pro- 
duced by a one-dimensional model containing 2 Mq of He and 
0.2 Mq of hydrogen. H and He are mixed and distributed out to 
~ 10000 km s _1 . A clumping factor of the He layer of £ = 200 
was used. As expected the [O i] AA 6300, 6363 doublet is single 
peaked. 



radioactive energy deposition as described by [Hou96] even 
if it was a SN lib and that some other mechanism is needed. 
To avoid a double-peaked profile of the [O i] AA 6300, 6363 
doublet, the hydrogen must be concentrated below 11000 
km s _1 . This is in perfect agreement with observations of 
the early-time Ha absorption minimum (see Table 3). 

SN 2007Y has the highest Ha to [O i] AA 6300, 6360 
doublet flux ratio (and also total flux) of all five SNe (see 
Figure [6} . Unless a serious amount of hydrogen has been 
accreted from a thick circumstellar wind, SN 2007Y must 
have had a significant fraction of H in its outer layers at the 
time of explosion. In addition, Ha is observed at ~ 15000 
km s during the first days after explosion. 

Therefore, SN 2007Y is most likely a SN of Type lib 
similar to SN 1993J and SN 2008ax. Then it is important 
to understand why the early time Ha absorption was weak 
i|Stritzinger et al.l 120091 ). A simple explanation may be that 
the hydrogen of SN 2007Y is more mixed with helium than 
the hydrogen of SN 1993J and 2008ax. In this case hydro- 
gen could be ionised more strongly at early times, which 
means that there is hydrogen but no H I. Thorough mixing 
of hydrogen and helium could also explain the strong, low- 
velocity Ha emission observed at late epochs (see Section 
[5]l . This is discussed in more detail in Section [6] 



5 AN ALTERNATIVE TO SHOCK 
INTERACTION 

To our knowledge there are only three stripped CC-SNe 
(SNe 1993J, 2007Y, 2008ax) which show strong, box- 
shaped Ha emission in their nebular phase . It is com- 
monl y assumed (e.g. Filippenko et al.l Il994l ; IPatat et al] 
1 19951 ; iHouck fc Franssonl ll996J ) that the box-shaped Ha 



emission of SN 1993J is caused by shock interaction. 

[Hou96] have shown that radioactive energy deposition 
is too weak to produce strong Ha emission at late epochs. 
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Figure 8. Optical spectrum of SN 2007Y at 288 days after ex- 
plosion (black line). The synthetic flux (red dashed line) was pro- 
duced by a one-dimensional model containing 1.5 Mg of He and 
0.1 Mq of hydrogen. H and He are mixed and distributed out to 
~ 10500 km s — 1 . A clumping factor of the He layer of £ = 100 
was used. As expected the [O i] AA 6300, 6363 doublet is single 
peaked. 
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Figure 9. Optical spectrum of SN 2008ax at 359 days after ex- 
plosion (black line). The synthetic flux (red, dashed line) was pro- 
duced using a one-dimensional model. A clumping factor of the 
He layer of £ ~ 100 was used. While the Ha emission is produced 
by tiny fractions of H in the He layer below 12000 km s , the 
double-peaked profile of the [O i] AA 6300, 6363 doublet is again 
produced by Ho absorption between 12000 and 12500 km s 
(also see Figure [3} • 



Their analytical estimates do not depend on the SN 1993J 
model explicitly and should be valid for any SN lib with 
moderate amounts of H (~ 0.1 Mg). One may therefore 
conclude that the same mechanism operating in SN 1993J 
is also at work in SNe 2007Y and 2008ax. However, for 
SN 2008ax there are several problems connected to the 
shock-interaction interpretation of its Ha nebular emission 
[Ta ulO]. In SN 2007Y shock int eraction is not expected at 
all (jChevalier fc Soderberdl2010r ). 

In this section we present an alternative mechanism for 
powering late nebular Ha emission. This mechanism may ex- 



plain how late-time Ha emission can be powered by radioac- 
tive decay without any need for shock interaction, solving 
the problems mentioned above. This mechanism is simple 
but needs special conditions to operate, which could explain 
why strong nebular Ha emission is rare in stripped CC-SNe 
(this may also be explained by shock interaction). 

We assume that a small fraction of hydrogen (of the 
order of 1%) is mixed into a helium-dominated layer. The 
degree of H ionisation will lie between 1 and almost 100% 
depending on the ionisation state of helium, since He re- 
combination can ionise hydrogen strongly (see Appendix A). 
There is a sharp transition between complete and almost no 
ionisation of hydrogen, if the hydrogen mass is comparable 
to the mass of He II. In such a situation, when the degree 
of ionisation of helium decreases with time, as it does nat- 
urally, large amounts of neutral hydrogen can be produced 
in regions where no H I had been present at earlier epochs 
(this is explained in more detail in Appendix A). Therefore, 
the H I fraction in the He dominated layer can increase from 
approximately zero to, say, 2% of the total mass within sev- 
eral hundred days. Up to this epoch an observer would not 
even know about hydrogen in these regions since H II can- 
not be observed, except by extremely weak recombination 
radiation. 

It is important to note that the energy deposited in the 
He-rich regions (at least in SN 2008ax, but probably also in 
SNe 1993J and 2007Y) is ~ 10 times larger than the energy 
emitted in Ha. This means that if only 10% of the deposited 
energy was transformed into Ha, shock interaction would 
not be needed to explain the late-time Ha emission. 

Under typical nebular conditions ~ 70% of the energy 
deposited in a He rich layer by radioactive decay is stored 
in thermal electrons. It is therefore clear that emission from 
this region is dominated by thermal electron collisional ex- 
citation and not by recombination. If heavy elements (C, 
O, Na, etc.) are mixed into the He layer in sufficient frac- 
tions, they can be excited by thermal electrons and cool the 
gas efficiently. In this case Ha is predominantly produced 
by recombination, which is much too weak to explain the 
late time Ha flux [Hou96; they neglected thermal electron 
excitation of hydrogen]. Since this scenario would produce 
high- velocity optical lines (e.g. oxygen or sodium), which are 
not observed, helium cannot be mixed with heavy elements 
strongly (at least in SN 2008ax; see Section [3}. 

If the gas consists of hydrogen and helium only, or if 
heavier elements are singly ionised to almost 100%, which 
would be the case if the fraction of these elements was lower 
than the fraction of He II, thermal electrons cannot cool 
effectively until they reach high enough temperatures to ex- 
cite hydrogen (~ 10000K ; see Appendix A) . In this scenario 
99% of the thermal electron energy is radiated away by Lya 
emission and therefore lost to the UV. About 1% is con- 
verted into Ha, which is too small by a factor of ~ 10 to 
explain the Ha observations. 

At this point clumping could become important. We use 
the symbol £ for the clumping factor, which is defined as the 
inverse of the filling factor [see iMaurer fc Mazzalil (|2010t ) for 
more details]. A clumping factor £ means that the density 
is increased by a factor of £ locally, while the global density 
remains constant. 

A helium layer containing small amounts of hydrogen 
can emit much more than 1% of its thermal electron energy 
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in Ha if clumping is strong, since clumping increases the 
excitation of H I to n = 2 & 3, and the self-absorption optical 
depth of the H I ground state transitions. 

Summarising, this means that, if hydrogen and helium 
are mixed in suitable fractions and clumped strongly, ra- 
dioactive energy deposition can power Ha completely with- 
out any additional source of energy. We present some esti- 
mates here, which support this idea. 

Since the excitation potential of hydrogen is ~ 10 eV 
lower than that of helium, the thermal electron col- 
lisional excitation rates of hydrogen are approximately 
exp(10eV/kT) ~ 10 5 " 3 times (at 10000 - 20000K) higher 
than those of helium. Unless there are ~ 10 5-3 times more 
He I than H I atoms ~ 70% of the total energy deposited in 
helium will be radiated away by hydrogen and He IR lines 
(the other 30% will go into non-thermal electron excitation 
and recombination radiation, mainly UV). 

It is clear that, as long as H is ionised completely (3> 
99%), Ha emission and scattering are extremely weak. There 
is no Ha scattering since there is no Lya self-absorption 
if there is no H I, which means that the H I n = 2 level 
is depopulated strongly. Moreover, there is almost no Ha 
emission, since only ~ 20% of the deposited energy go into 
ionisation. This energy is mainly emitted by He I recom- 
bination radiation and Lya, which means that the fraction 
of deposited energy emitted by Ha is less than one percent. 
Thermal electrons cool predominantly by exciting He I, since 
the temperature increases to more than 15000K. 

We derive a rough estimate for the dependence of Ha 
emission on clumping. We consider effective recombination 
into the hydrogen levels n = 2 and 3, thermal electron (de-) 
excitation between the ground state and the levels n = 2 
and n = 3, and radiative transitions from H I n = 3 to n = 
2 and from n = 2 to n = 1. We neglect radiative transitions 
from Hln = 3ton = l, since this rate is small owing to self- 
absorption. We also neglect thermal excitation from H I to n 
= 2 to n = 3, which may contribute to the population of the 
H I n = 3 level but is less important than direct excitation 
from the ground state. We only consider the situation where 
enough H I is present to cool the gas efficiently ( H[I1 > 
0.1%), which means that the H I ground state transitions 
are highly self-absorbed. 

To calculate the formation of Ha the population of the 
n = 2 state of H I has to be estimated 
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the n = i and j states of H I (e.g. Callaway 1994), I is the total 
ionisation rate (non-thermal electrons and photo-ionisation; 
see Appendix A), f n are the fractions of recombining elec- 
trons, cascading into the H I n states and 
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is the Sobolev optical depth of the H I ground state, with 
ALya ~ 1216 A the wavelength of Lya, g n the statistical 
weights of the n states of H I and A21 the radiative rate of 
Lya. Depending on the number density of the H I ground 



state the n = 2 state can be thermally populated (C21C, 2> 
■^t) at moderate clumping factors (£ ~ 10) already. The 
Ha optical depth is given by 
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where \n a ~ 6563 A is the wavelength of Ha and A32 is 
the radiative rate of Ha. th q increases with clumping at 
low clumping factors, but then saturates as soon as the de- 
population of the H 1 n = 2 level is dominated by collisional 
de-excitation. In rough approximation the ratio of the Ha 
to the Lya luminosity is given by 

Ha -Ehq n H[i],3 A' 32 
Lya Eh ya n H [i],2 A' 21 
E Ha gigCj+J/g] 



El^ya Cl2C[+-f/2 
C21T21C 
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A21 JV ' A 32 [l -exp(-THaC)] 

~0.03-exp(-^), C~l,Cij>/ 
~ F(T, n ,n H ,n H [ii],C), C>1 

where A'^ denotes the self-absorbed radiative transition 
rates between the n = i and j states of H I. 

Clumping can increase the relative strength of Ha sig- 
nificantly. High clumping factors (£ 2> 1) increase the Ha 
optical depth less effectively than they do the Ha luminos- 
ity, while moderate clumping factors increase both. An es- 
timate for the time-dependent density n H [i] is available (see 
Appendix A). 

A numerical computation of the processes described 
above can be performed using our nebular code. In Figures 
El E] and H] we show models of SN 1993J (363 days), 2007Y 
(288 days) and SN 2008ax (359 days) using the mechanism 
described above (radioactive decay energy; mixing of H and 
He; strong clumping). Since the late Ha emission of SNe 
200 lig and 2003bg is much weaker it is clear that the late 
Ha emission of these SNe can be reproduced i n a similar way. 
For SN 1993J we use the He mass estimate of lWooslev et al.l 
l|l994h . For SN 2007Y we assume that the He mass is a bit 
smaller than in SN 1993 J. For SN 2008ax the He mass is 
estimated in Section [3] For SN 1993 J we use a helium mass 
of 2 Mq, a hydrogen mass of 0.2 Mq and a clumping factor 
C ~ 200, for SN 2007Y we use a helium mass of 1.5 M , a 
hydrogen mass of 0.1 Mq and a clumping factor C, ~ 100 and 
for SN 2008ax we use a helium mass of 2 Mq, a hydrogen 
mass of 0.2 Mq and a clumping factor of £ ~ 100. Most of 
the hydrogen is concentrated in a thin shell at high veloc- 
ities, while tiny fractions are mixed into the lower velocity 
helium. It is difficult to decide whether clumping, the helium 
or the hydrogen mass should be increased to obtain a repro- 
duction of the spectra. Several combinations are possible. 
Additionally, the outer region could be influenced by shock 
interaction. Therefore it is not possible to derive a hydro- 
gen or helium mass estimate from this procedure. It seems 
however clear that small amounts of hydrogen are sufficient 
to reproduce late Ha emission (which is usually not possible 
[Hou96]), if there is strong clumping (£ ~ 100). 
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It is important to note that we could not reproduce 
the temporal evolution of the Ha line using this mechanism 
with our nebular code. We can find models that reproduce 
the observations at any epoch, but no model describing the 
observations at all epochs consistently. This is discussed in 
Section |S1 



6 DISCUSSION 

6.1 Late Ha emission 

We have investigated the formation of Ha in the nebular 
phase of five SNe lib (re- classifying SN 2007Y as SN lib) . We 
find that radioactive energy from 56 Ni decay is not sufficient 
to power the Ha emission observed at late nebular epochs 
(> 200 days), as long as the energy deposited in the He layer 
cannot be tapped. This finding is consistent with the work 
of [Hou96] for SN 1993J, but raises the question of which 
mechanism causes the observed Ha flux. 

In Section [5] we have shown that a combination of mix- 
ing and strong clumping could solve this problem. It is not 
clear whethe r such high clumping factor s are realistic (how- 
ever see e.g. iKozma fc Franssonl (1998)). Furthermore, we 
could not reproduce the temporal evolution of Ha using this 
mechanism. This may be because of several reasons. First, 
there is a degeneracy between the hydrogen and the helium 
mass and clumping and it is not clear which combination has 
to be chosen. Secondly, as can be seen from Equationf5J the 
strength of Ha is very sensitive to the electron density and 
temperature at clumping values of the order 100. In addi- 
tion, n H [i,n] is extremely sensitive to mixing and clumping. 
The composition has to be known at the one percent level or 
better. Our nebular code is expected to be less accurate. We 
do not treat UV-radiation transport explicitly, the ionisa- 
tion of hydrogen by helium is treated in approximation only 
and in general there are uncertainties owing to the atomic 
data and owing to the incomplete description of the physical 
scenario. Asymmetries of the ejecta can further complicate 
the problem. 

Since the radioactive scenario appears to have problems, 
it may be worth considering the most common interpreta- 
tion, i.e. that late SN lib Ha emission is the result of shock 
interaction. An extensive discussion of several shock interac- 
tion scenarios for SN 2008ax can be found in [TaulO]. How- 
ever, as pointed out by [TaulO], there are several problems 
within all these scenarios. Shock interaction has difficulties 
explaining why Ha emission is observed at low velocities in 
SN 2007Y and SN 2008ax (see [TaulO]). Moreover, there 
are contradictions between X-ray and Ha observations, if 
both are interpreted as shock interaction (see SNe 2003bg, 
2007Y). These problems may arise from an incomplete un- 
derstanding of the influence of the shock radiation on the 
ejecta or from an inaccurate treatment of the physical pa- 
rameters determining the shock properties (e.g. the SN outer 
density structure). 

The shock interaction model can explain why strong Ha 
emission at late times is rare in SN lib and has the advantage 
that it may explain a flattening of the late Ha light curve, 
which is observed in SN 1993J and possibly in SN 2008ax. 
However, this has never been shown quantitatively. Nor has 
it been shown that shock interaction can reproduce late Ha 
emission in detail in any SN lib at all. 



The mechanism presented in Section [5] can explain the 
formation of low velocity Ha in SNe 2007Y and 2008ax, 
which is observed. It can further explain the absence of 
high velocity Ha emission, again consistent with the ob- 
servations, since in this scenario the Ha emission traces the 
He distribution. It can explain why SN 2007Y, has strong 
Ha emission, although IChevalier fc Soderberd (J2010J) claim 
that there is extremely weak shock interaction in SN 2007Y. 
It can ex plain why SN 2003bg has no strong Ha emission, 
although Chevalier fc Soderberd 1)2010] ) claim that there is 
strong shock interaction in SN 2003bg. Since the model 
needs some fine tuning it can explain why strong late time 
Ha is rare in SNe lib. 

There may well be some other mechanism (not clump- 
ing; still mixing of H and He), which emits in Ha energy 
absorbed by the He layer. If it was possible to use ~ 10% 
of the energy of a massive He shell to excite Ha without 
producing any other strong optical lines, the late-time Ha 
emission could be explained by radioactive energy deposi- 
tion alone. 



6.2 Nebular line profiles of SNe lib 

Usually, double-peaked line profiles are interpreted as 
toroidal ejecta distributions. While we agree with this in- 
terpretation in general (e.g. for SN Ib/c; see also Section 
[3} we think that the situation may be different in SNe lib. 
We have shown that, taking the 40 day (after explosion) Ha 
absorption minimum velocity as lower boundary of the bulk 
of hydrogen (as it is usually done) we can explain the peak 
profile of the [O i] AA 6300, 6363 doublet for all five SNe of 
our sample consistently by Ha absorption. 

Since Ha absorption causes a split of the [O i] A A 6300, 
6363 doublet if it is located around 12000 km s , and since 
the position of Ha influences the position of this split, we 
can infer the radial distribution of H I from fitting the [O i] 
AA 6300, 6363 doublet profile. For all five SNe of our sample 
the Ha velocity measured with this method agrees very well 
with the early-time Ha minimum velocity (see Table 3). It 
seems unlikely that this is coincidence. 

The outer Ha region could be optically thick even at 
late epochs (e.g. excited by shock interaction), which could 
explain why the profile of the [O i] A A 6300, 6363 doublet 
does not change its shape significantly at late epochs. For 
example, the reverse shock of SN 2008ax could be located 
around 12000 km s _1 at 359 days after the explosion [TaulO]. 



6.3 SN 2008ax 

We have derived a three-dimensional model of SN 2008ax. 
Two kinds of asymmetries may be observed. First the He IR 
lines are double-peaked, which is most easily explained by 
a torus shaped He distribution. Secondly, the blue sides of 
the He lines appear to be stronger than the red sides. 

We find a total ejecta mass of ~ 3 Mq , containing about 
0.1 Mq of 56 Ni and expanding with a kinetic energy of ~ 
10 51 ergs, which is consistent with findings from light curve 
and early-time Ha absorption modelling (Table 2; also see 
[TaulO]). We find a helium mass of ~ 2 Mq and derive the 
abundances of heavier elements such as carbon or oxygen 
(Table 1). 
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While most of the heavy element mass is confined to ve- 
locities below 4500 km s _1 , helium is predominantly found 
between 4500 and probably 12000 km s _1 . The outer bound- 
ary is uncertain. Some helium is present below 4500 km s _1 . 
We find that a separation of helium and heavy elements, at 
least between 2000 and 4500 km s _1 , improves the repro- 
duction of the He IR lines. 

Since our nebular code does not reproduce the contin- 
uum flux observed around 20000 A, we added some artifi- 
cial background to allow He I 20587 A line scattering. There 
seems to be no He I 20587 A line scattering above ~ 10000 
km s _1 , while the He I 10830 A line would not be in conflict 
with higher helium velocities. The observations of the He I 
20587 A line place an upper limit on the He mass between 
9500 and 12000 km s" 1 of ~ 0.05 M . ft may well be that 
continuum destruction (e.g. by hydrogen), which could re- 
duce the He I 20587 A optical depth, is important. In this 
case the He mass in this region could be larger. 

We can compare our results to SN lib m odels 
of iNomoto et all (119931 ) and IWooslev et all (J1994T) , The 
iNomoto et alj ( 19931 ) models [also see Hou96] predict a 
strong confinement of 56 Ni to low velocities (< 2000 km 
s _1 ). Other heavy elements are distributed in layers above 
the radioactive core. Most of the hydrogen is confined to 
a thin layer above the He layer. This He layer contains 
some traces of hydrogen a nd heavier elements. The mod- 
els of IWooslev et al.l |l994j) predict a broader distribution 
of 56 Ni in velocity space. Heavy elements dominate at low 
velocities but are also mixed into the He layer. Helium and 
hydrogen are mixed in the outer regions. 

We find that 56 Ni is probably not as centrally confined 
as predicted by the INomoto et al.l (1 19931 ) models. On the 
other hand carbon seems to be confined to a thin shell on top 
of the iron core in agreement with these models. Since our 
core is deformed to an axis ratio of ~ 3:2 (see Figure[T}, one 
probably cannot expect perfect agreement with the spherical 
symmetric models. The distribution of heavy elements seems 
rather consistent with the predictions in the core as well as 
in the He layers. Our results for hydrogen are inaccurate and 
therefore a comparison is difficult. 

We compare the masses of He, C, O, Ca and 56 N i of ou r 
SN 2008ax model to the 13C model oflWoosley et all jl994h . 
which has been used bv lTsvetkov et al.l ( 20091 ) to reproduce 
the light curve of SN 2008ax. We find that our estimates of 
the He, C, Ca and 56 Ni mass agree at the 10% level. The 
estimate of the oxygen mass is larger by a factor of ~ 2 than 
predicted by the 13C model. 

Thanks to its late time IR observations, SN 2008ax is 
the first SN lib where the He density field can be deter- 
mined from nebular observations. Other methods, such as 
light-curve modelling, provide poor information about ele- 
ment abundances or asphericities. Early-time modelling can 
determine the properties of the outer layers more accurately, 
but no information about the central region can be obtained. 
Therefore nebular IR observations (10000 to 22000 A) of SNe 
lb and lib are highly desirable. 



7 SUMMARY & CONCLUSION 

We have derived estimates related to hydrogen and helium 
line formation in the nebular phase of SNe lib. We have 



shown that hydrogen can be highly ionised if it is mixed with 
helium in suitable ratios. This has probably been known 
before. However, it could explain why Hq can sometimes be 
observed in the nebular phase but not at early times. This 
finding is of special interest for SN 2007Y. 

From our analysis it seems likely that SN 2007Y is a 
SN lib. The lack of strong Ha absorption around 50 days 
after explosion may well be explained by strong mixing of 
H and He. Such mixing could also explain the formation 
of strong late time Ha. Since we cannot determine the He 
distribution of SN 2007Y, this remains a speculation. 

We have shown that radioactive energy deposition is 
insufficient to produce the late-time Ha emission in all SNe 
lib of our sample, as long as the energy of the helium layer 
cannot be tapped by Ha. 

Shock interaction may be observed in all SNe of our 
sample. There seem to be two types of late Ha emission 
scenarios. Strong, box-shaped (SNe 1993J, 2007Y, 2008ax) 
and weak (SNe 2001ig, 2003bg) Ha. The simplest expla- 
nation would seem to be strong and weak shock interac- 
ti on, respectively. Howeve r , this is in conflict with findings 
of IChevalier fc Soderbergl (J201Ch . at least for SNe 2003bg 
and 2007Y. It also seems difficult to explain the low velocity 
observations in Ha by shock interaction [TaulO]. 

We presented an alternative mechanism to shock in- 
teraction, explaining late time Ha emission by radioactive 
energy deposition. The right combination of mixing and 
clumping of hydrogen and helium has been shown to be able 
to reproduce the Ha observations of all SNe of our sample at 
least up to 350 days after explosion. However, there are also 
some problems with this interpretation of late Ha emission. 
Summarising, it is not clear how late Ha is formed in SNe 
lib. Therefore, we were not able to derive an estimate of the 
hydrogen mass for any SN of our sample. 

We have shown that most likely the profile of the [O 
i] AA 6300, 6363 doublet is influenced by Ha absorption 
strongly. Hydrogen concentrations above ~ 11000 km s _1 
can cause double-peaked oxygen profiles, while slower hy- 
drogen can not. We have shown that observations of early- 
time Ha absorption minima and the corresponding oxygen 
line profiles are perfectly consistent with this interpretation 
for all five SNe of our sample. 

This scenario also explains why other lines, like Ca II or 
[O i] 5577 A, are single peaked in all SNe of our sample. An 
exception are the He lines of SN 2008ax, which show signs of 
asymmetry. It seems likely that at least the inner part of SN 
2008ax is asymmetric. We have obtained a three-dimensional 
model of the SN 2008ax envelope. For the first time the he- 
lium mass of a SN lib has been determined from nebular 
modelling. We have obtained estimates of the total mass and 
kinetic energy of SN 2008ax in excellent agreement with re- 
sults from light-curve modelling. Chemical abundances have 
been derived. The nebular model of SN 2008ax provides the 
opportunity to compare observations and theoretical SN lib 
models with unpreceden ted richness of detai l. 

As already stated bv lHouck fc Franssonl |l996|), nebular 
IR observations (10000 to 22000 A) of core-collapse SNe 
would be extremely useful for deriving properties of CC-SNe 
such as helium abundance, total mass and kinetic energy. It 
is important that both He IR lines (10830 and 20587 A) 
are observed simultaneously. Such observations should be 
obtained for any SN lb or lib. 
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APPENDIX A 

Hydrogen ionisation in a He-dominated layer 

We derive some estimates for the ionisation of hydrogen by 
non-thermal electrons and by the UV-radiation emitted by 
helium. We assume that H is mixed with He but not with 
other elements. Therefore the treatment presented is only 
valid if heavy elements are present in smaller fractions than 
hydrogen and helium. However, other elements with ionisa- 
tion potentials below ~ 20 eV would be affected in a way 
similar to hydrogen. 

About 40% of all helium recombinations go directly to 
the ground state and cause another ionisation of helium and 
possibly of hydrogen. This recycling process increases the 



© 0000 RAS, MNRAS 000, 000-000 



14 Maurer et al. 





(J 

-1 


^v 


I 


-2 




- 


X 


-3 




- 


o 
O 


-4 




- 




-5 




^^____^ 




-R 


~~ 



3 4 5 6 7 

He[ll]/He [%] 



'0 




1000 
Epoch [d] 



/000 



Figure 10. The logarithmic fraction of neutral hydrogen —3111 
(full black line; Equation IIOI I as a function of the He II fraction 
H< 4 n l [%] for a helium dominated layer containing 2% of hydro- 
gen, assuming / = 0.5 and §H = 0.1. The He density was set 
to 10 8 g cm -3 . There is a steep decrease of H I as soon as the 
fraction of He II increases above -j%. Changing the degree of ion- 
isation of He by 1% only (from 4.5 to 3.5% in this example), the 
fraction of H I increases by a factor of 10 4 . 
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Figure 11. Toy model of a SN lib. The time-dependent fraction 

H< 4 n l [%] (full black line; Equation I15H is shown for a model 

described in the text (M He = 3 Mq, M Ni = 0.1 M Q , / = 0.5, T 



15000K). The fraction 



"H[I] 



[%] (dashed red line; Equation 



»H+»Ho 

I10|l is shown for hydrogen masses of 0.01 to 0.1 Mq (bottom to 
top). There is no clumping of the H and He layer. The calculation 
starts at 100 days since at earlier epochs photospheric radiation 
may be important for ionising hydrogen and helium and our esti- 
mates become invalid. Although the model is highly simplified it 
demonstrates that depending on the ratio of hydrogen and helium 
the H I fraction can be approximately zero for the first 100 days 
and can increase rapidly at later epochs. This means that no Ha 
would be observed at early epochs, while strong Ha emission or 
scattering is possible later. 



Figure 12. Toy model of a SN lib. The time-dependent fraction 
— il— L [%] (full black line; Equation 1 15D is shown for a model 
described in the text (M Hc = 3 Mq, M Ni = 0.1 M©, / = 0.5, T 



15000K). The fraction 



"H[I] 



[%] (dashed red line; Equation 



»H+»Ho 

110ft is shown for hydrogen masses of 0.001 to 0.01 Mq (bottom 
to top). The clumping factor of the H and He layer is set to 100, 
which reduces the degree of ionisation strongly. The calculation 
starts at 100 days since at earlier epochs photospheric radiation 
may be important for ionising hydrogen and helium and our esti- 
mates become invalid. Although the model is highly simplified it 
demonstrates that depending on the ratio of hydrogen and helium 
the H I fraction can be approximately zero for the first 100 days 
and can increase rapidly at later epochs. This means that no Ha 
would be observed at early epochs, while strong Ha emission or 
scattering is possible later. 



total ionisation rate but is not important for the rest of the 
appendix. 

Electrons recombining to an excited state reach the 
ground state by two-photon emission (2PE) of the 2s( 1 S) 
state during the nebular phase. At early times, when the 2p 
levels are strongly excited by radiation, most electrons reach 
the ground state via the 2p levels. 

In the case of 2PE approximately 30% of the ground 
state transition radiation can ionise hydrogen (the two pho- 
tons are created with energies between and 20.6 eV; the 
chance of produci ng a photon with E > 13.6 eV is ~ 30%; 
iDrake et al.lll969l ) and in the case of 2p transitions 100%. 
Therefore, we assume that a fraction / (0.3 — 1.0) of all 
He recombinations into excited states can ionise a hydrogen 
atom. 

The total ionisation rate of hydrogen is then determined 
by non-thermal electron ionisation and by radiative ioni- 
sation from He recombination radiation. The non-thermal 
electron ionis ation rates of n eutral hydrogen and helium are 
given by (e.g. lAxelrod|[l980l ') 



y = K 



H[I],Hc[I] 



^Dcp 



A r TotW / H[I],Ho[I] 



(5) 



where Nxot is the total number of atoms and WH[i],He[i] is 
the work per ion of hydrogen and helium, which depends on 
the ratio of electrons to atoms and on the absolute atomic 
density (weakly). 

For simplicity we assume W^nm = M'Hcfii in this section, 
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which will cause an error of ~ 30% on the ionisation rates 
of H and He in the worst case. 

During the fi rst few hundre d days the SN is in ionisation 
equilibrium (e.g. lAxelrod 1980) and the ionisation balance 
of hydrogen can be estimated as 

Y(n He [i]D Hm f + n H [i]) = #h[i] («-hc[ii] + WH[n])C n H[n] (6) 

where -D H [i]/ is the fraction of He recombination radiation 
ionising hydrogen with 



A 



1 — cxp(— (TnARrijin 



(7) 



^H[i] = J. - exp^-CTH^-n,'«H[i]J 

where <th is the ionisation cross-section of H at ~ 20 eV 
and AR is some characteristic width of the He shell. The 
ionisation balance of Helium can be written as 

Ynn e [I] = Rne[I](nHe[lI] + ™H[II] )C M Hc[II] (8) 

Setting R = -Rh[i] = ^Hc[i] (which is a good approximation) 
one obtains 

"•Hem 



-Dh[I]«Hc[I]/ + «H[I] 



(n H -n H [i]) 



"•Ho[II] 

The exact solution of Equation [9] is given by 



(9) 



._!„./ n H e[I]«He[II] ^ , [ ™Hc[I] £h ™H - /n H e[II] ) 

H ("h - /n Hc [ii]) 

"He 

(10) 

where W{x) = X^^Li , — x n is the Lambert W function. 

In approximation, the deposition fraction -Dmi] is given by 



£>H[I] = 



1 



<T H Ai?n H[I] > 1 
<THAi?n H [i] < 1 



<T H Ai?nH[i] = {h«h[i] 

The ionisation cross section of hydrogen is ~ 2 -10 -18 cm 2 at 
20 eV and the SN radius is of the order of 10 16 cm during the 
first few hundred days, which means that £h will be of order 
0.01 — 0.1 cm 3 . The number density of neutral hydrogen is 
then given by 






"Hs[I] ( 
"Ho \ 



1-/ 

Ho[I] 



'Hc[II| 



Ch^ H [I] > 1 

(fn He {u]^n) CHn H [i] < 1 



l_ nHc+»lHo[I]"Ho[II]/?H 

The degree of ionisation of hydrogen increases rapidly as 
soon as the hydrogen number density becomes similar to 
the number density of He n times / (see Figure [T0|l . 

A factor of 10 or less in the hydrogen abundance can 
make a difference of several orders in magnitude in the de- 
gree of ionisation. Also, a small change in the ionisation 
balance of He can have serious influence on the ionisation 
balance of H. A simple example is given by a helium layer 
containing 2% hydrogen and %^% He II. Hydrogen is ionised 
to 3> 99% then. If the amount of He II decreases to lj-% (a 
decrease of the ionisation fraction is expected at late epochs; 
see Figures [TTl and [T2"Tl the H II fraction decreases to less than 
10%. Without changing the physical conditions sig- 
nificantly, the fraction of H I can change by several 
orders in magnitude. 



Toy model of a SN lib 

To obtain some rough estimate of the electron tempera- 
ture in a He-dominated layer containing some small fraction 



of hydrogen, we calculate the temperature-dependent ratio 
of collisional excitation and recombination 1Z. The thermal 
electron excitation coefficient from n = 2 to n = 1 is given 
by 



C12 ~ 8.6 x 10~ 
= Con e 
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(11) 



where f^i is the effective collision strength of the hydr ogen 
n = 1 to n = 2 transitions given bv lScholz et al.l (| 19901 ) and 
/2,<?i are defined in Section [5] This gives 
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(12) 



While Co increases with temperature 7? decreases and -^j- 
is approximately one at ~ 11000K. Since hydrogen must be 
excited by thermal collisions effectively in order to produce 
a significant contribution to the total luminosity, 1Z must be 
larger than 2> 1. Therefore the temperature must be higher 
than ~ 10000K, depending on the H I fraction. 

At temperatures of ~ 20000K (which are reached for low 
fractions — — < 0.1; this implies nnm ^ ^Hefii] ^ n Hc[i]) 
thermal excitation of He becomes important, the above es- 
timate becomes invalid and the relative importance of H I 
emission decreases as compared to He I. Therefore, the tem- 
perature should be ~ 10000 — 15000 K in order to allow 
effective H I emission. 

We can use this temperature estimate to calculate the 
temporal evolution of He n and therefore H I for a toy model 
of the SN ejecta. Of course this temperature estimate is not 
exact, but it is sufficient to demonstrate how the absolute 



fraction 



•*H[I] 



in a helium dominated layer can evolve with 



»H+"Hc 

time (see Figures HTl and IT2T) . 

The deposited energy can be estimated analytically if 
the He density is known. For simplicity we assume that X 
solar masses of He and Y solar masses of 56 Ni are distributed 
within a sphere of velocity v homogeneously (of course this is 
in conflict with our assumption that He is mixed with small 
fractions of heavy elements only, but this is a toy model 
only and He and 56 Ni could be separated on small scales). 
At nebular epochs (say, > 100 days), when most of 56 Ni has 
decayed to 56 Co, the deposited luminosity is given by 

L Dcp {t) ~ 1.3 xlO 43 exp (-3^)M5 6Ni [M ]/ Dcp ergs s" 1 

(13) 
where /d cp is the deposition function, which can be calcu- 
lated analytically for a homogeneous sphere (e.g. lAxelrod! 
Il980f ). For three solar masses of helium and a 5 Ni mass 
of ~ 0.1 Mq distributed homogeneously in a sphere of 
10000 km s _1 the luminosity deposited around 350 days is 
~ 10 39,5 ergs s" 1 while the Ha luminosity of SN 2008ax at 
this epoch is ~ 10 38,5 ergs s _1 only [TaulO]. 

Since the ionisation potential of He is high and since 
we assumed that it is the predominant element, helium is 
mainly ionised by non-thermal electrons 

TZH e YHcn Ucm = -RhcC(?1Hc[II] + "H[II])"Ho[II] (14) 

where 7£hc ~ 1-5 is the recycling fraction of helium. This 
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gives 

™Hc[II] ^HcC n H[II] + T^-Hc^Hc 

flHe 2_RhoC«-Ho 



(i?HcC n H[II] + 'R-UcYhc)' 



(15) 



We compute Equations \W\ and [15] for the toy model de- 
scribed above. We assume that 50% of all He recombinations 
into excited states can ionise a hydrogen atom (/ = 0.5) 
and vary the hydrogen mass between 0.01 and 0.1 Mq. The 
temporal evolution of the H I fraction is shown in Figure [TT1 
Depending on the ratio of hydrogen and helium the fractions 
of H I can increase dramatically between 100 and 1000 days. 
This means that Ha scattering or emission can appear at 
late epochs, without any hydrogen being detected at earlier 
times. In Equation \W\ and [15] we assumed ionisation equi- 
librium, which may become invalid at epochs of ~ 500 days 
and later. Therefore, we underestimate the He [II] fraction 
at very late epochs. However, qualitatively the effect is the 
same, since H I ionisation is related to He II recombination 
even at 1000 days and later. 
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